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Memory functions, which enter into the equations of motion for time correlation functions, are 
constructed from neutron scattering, infrared absorption and light scattering data involving single 
particle motions in liquids. The qualitative features of these memory functions are related to the 
shape of the corresponding time correlation functions. It is found that a negative portion to the memory 
function is indicative of a rapid loss of correlation in time while strong temporal correlations imply a 
memory function which does not go negative. The mathematical structure of a memory function is 
examined for the case of the ideal gas by expanding and evaluating the projection operator representa- 
tion of the function. The resulting expression has a rich mathematical structure and can be expressed 
in a closed form only for its Laplace transform. 
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1. Introduction 

The time correlation of molecular quantities has 
proven to be a useful way of describing molecular 
motions in liquids. These time correlation functions, 
which are obtained as Fourier transforms of data, 
can be investigated using recently developed techniques 
of time dependent statistical mechanics [1]. Exact 
equations of motion for the time correlation functions 
have been derived which involve memory functions, 
quantities which characterize the history of the fluid 
[1, 2]. While time correlation functions are now widely 
used to represent molecular motions [3], only a few 
studies of memory functions have appeared [4-9]. 
In this paper we examine some properties of memory 
functions as derived from data and from model cor- 
relation functions. Our results, combined with other 
results obtained from computer simulations, enable us 
to make some qualitative observations about the general 
features of these important quantities. 

A motivation for studying memory functions is the 
hope that these quantities will have simpler structure 
than the time correlation functions. Should this be 
true, it would be possible to model molecular motions in 
liquids in terms of a relatively simple memory function 
once the qualitative features of the memory functions 
are understood. Indeed, this has been done empiri- 
cally for the velocity time correlation functions 
obtained by computer simulation of the Lennard- 
Jones model fluid [10]. 

The mathematical structure of memory functions 
has been studied in the short time and long time limits 
[4-6]. We present here a calculation for all times of the 



memory function for the incoherent neutron scattering 
by an ideal gas. The mathematical structure of this 
quantity is found to involve even powers of the velocity 
of a particle in terms of deviations of the even powers 
from their average values. This structure is more com- 
plex than the structure of the time correlation function. 
Thus we have examples which suggest the existence of 
a rich underlying mathematical structure which yields 
relatively simple functional forms for the memory 
functions. 

It is plausible to expect that future studies will find 
a whole range of cases and that for as yet to be deter- 
mined physical situations, the memory functions are 
comparatively simple quantities. 

The paper is organized in the following way. First 
we introduce some notation and outline briefly the 
relation connecting a time correlation function and a 
memory function. In section 3 we construct memory 
functions from data and discuss the structural features 
of these quantities. The mathematical structure of 
memory functions is probed in section 4 where the 
memory function for the incoherent scattering from 
an ideal gas is calculated from the microscopic defini- 
tion of a memory function. The structure discovered in 
this analysis is found to have some unexpected features. 
The discussion of the qualitative features of memory 
functions in section 5 makes use of both our results and 
the results of other studies. 

2. Formal Preliminaries 

Molecular motions in liquids may be studied ex- 
perimentally in a number of ways. For example, the 
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translational motion of a single atom is probed by 
incoherent neutron scattering and the changes in 
orientation of a single molecule is probed by depolar- 
ized Rayleigh scattering. Infrared absorption and 
certain types of depolarized Raman scattering also 
probe orientational motions. In each case one can 
obtain from the data, a time correlation function 
which describes some aspect of the motion of a 
typical molecule in the fluid. 

The notation we shall employ for a time correlation 
function is \jj(t) with i//(0) = 1 where 



,M0 = <«*(0M0>/<M 2 >. 



(2.1) 



Here the angular brackets (...) indicate an 
ensemble average and the phase function 



u(t) = exp (-iLt)u(O) 



(2.2) 



evolves in time as determined by the Liouville operator, 
L, for the system. The phase functions associated with 
various experimental probes are described elsewhere 

[ii]. 

A formally exact equation of motion for i//(f) may 
be derived using the projection operator formalism of 
Zwanzig [1, 2]. The equation is 



#(*) 

dt 



-/: 



with 



ch// 
dt 



dt'K(t-t')x\s(t') 



= 0. 



(2.3) 



(2.4) 



The memory function K(t) is specified as 

K(t) = (u*(Q) exp [-(l-P)iLt]u(0)) (2.5) 
where the projection operator P has the property 

Pg=u(u*g) (2.6) 

for any phase function g and 

u = — iLu. (2.7) 

3. Memory Functions 

In this paper we examine memory functions ob- 
tained by recursively solving eq (2.3) for K(t) when 
i//(j) is known from experiment. The numerical pro- 
cedures are discussed in the appendix. We use in- 
coherent neutron scattering data for liquid argon 
[12] as an example of translational motion. Infrared 
absorption data and depolarized Raman scattering 
data for neopentane [13] and depolarized Rayleigh 
scattering data for benzene [14] are used as examples 
of orientational motion. 

The results are displayed in figures 1-4. There the 
correlation function, ^(t), and the normalized memory 
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FIGURE 1. Single particle translation motion for liquid 
argon at T = 85.2 K. 



The time correlation function (- 
= 2.0A _1 are pictured. 



and normalized memory function (-- 



for Q 




FIGURE 2. Reorientation of liquid benzene molecules at t = 
studied via depolarized Rayleigh scattering. 



59 °C 



The 1 — 2 component of reorientation motion is pictured. The time correlation function 
( ) and the normalized memory function ( ) are shown as functions of time. 

function, K(t)IK(0) , are shown as functions of time. 
These memory functions are characterized by a rapid 
initial decrease to negative values followed by an ap- 
proach to zero which is prompt when compared with 
the decay of the correlation function. This appears to 
be a fairly common situation as evidenced by its 
presence in the memory function for the incoherent 
neutron scattering from an ideal gas as shown in 
figure 5. It is not a universal feature as counter ex- 
amples may be constructed. 

The results for K(t) become erratic for long times. 
This is a result of the sensitivity to irregularities in 
ifj(t) of the numerical procedure used to obtain 
K(t). The uncertainty of K(t) in figure 1 is on the 
order of ±0.02 for times greater than 10~ 12 s. The pair 
of dashed lines in figure 3 indicate the range of values 
assumed by K(t). In figure 4 the solid horizontal lines 
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FIGURE 3. Reorientation of neopentane at 241 K studied by infrared 
absorbtion. 

The/= 1 component of reorientation motion is pictured. The solid line ( ) is the time 

correlation function and the dashed line ( ) is the normalized memory function. 
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FIGURE 5. The time correlation function ( ) and normalized 

memory function ( ) for incoherent neutron scattering from 

an ideal gas. 

The time is measured in units of q-I'lfSm. 




TIME 



FIGURE 4. Reorientation of neopentane at 239 K studied by Raman 
scattering from the 924 cm~ x band. 

The time correlation function ( ) and the normalized memory function ( ) for 

the 1 = 2 component of reorientation motion are displayed as functions of time. 



FIGURE 6. The model time correlation function ( ) and the 

corresponding normalized memory function ( ) of eq (3.1) 

are displayed. 

This shows that memory functions need not have negative regions. 
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indicate the range of values obtained for K{t) for times 
greater than 0.8 X 10~ 12 s. The data used to obtain 
figure 2 were carefully smoothed so the erratic behavior 
does not show up until longer times are reached and 
would not be visible with this scale. 

Figure 5 shows the memory function associated with 
a Gaussian time correlation function (ideal gas). 
Here the negative region is not as pronounced as in 
liquid argon but it is definitely there. 

A model time correlation function which does not 
have a negative portion to its memory function is 



i//(0 = exp [A - VA* + t 2 lr 2 ]. 



(3.1) 



The results for A = 1.0 and 7=0.5 are displayed in 
figure 6. The tail of the memory function is approxi- 
mately exponential. 

The existence of a negative portion in the memory 
function is not surprising. From eq (2.5) we observe 
that C(t), the negative of the second time derivative 
of ifj(t) 9 is not dissimilar to K(t) since 



C(t)=-d 2 ^(t)/dt 2 

= (a* exp [— iLt]u). 



(3.2) 



At least for short times, K(t) and C(t) should be 
fairly close in value. They are identical at t = 0. The 
function C(t) necessarily has negative portions since 



/; 



dtC(t) = 0. 



(3.3) 



This follows readily from the initial condition d\jj(t)l 
dt\ t = o = 0. By taking Laplace transforms we obtain 



H 



C(z)=\ dte~ zt C{t) 



(3.4) 



-z 2 i//U)+z 



so that 



/; 



dtC(t) = lim C(z)=0. 



(3.5) 



Let us now examine further the relationships between 
C(t) and K(t). The operator identity 



exp[— (1— P)iLt] = exp[ — iLt] 



+ (' dt'exp[-iL(t-t')]iPL exp[-i(l-P)Lt'] 
may be used to show that 



(3.6) 



K(t) = C(t) 



f'dt'K(t-t' 
Jo 



)d^(t')ldt' . 



(3.7) 



This equation may also be obtained directly from eq 
(2.3). In some of the cases considered in this article, 



the second term in eq (3.7) is numerically small com- 
pared to K(t) for short times. This suggests that 
K(t) ~C(t) might be a useful starting point for an 
iterative solution of eq (2.3). Unfortunately this is 
not the case as only the trivial solution ifj(t) = 1 exists 
whenK(j) is replaced by C(t). 

Another limitation of this way of examining K(t) 
only becomes apparent after the mathematical 
structure of K(t) is studied by evaluating eq (2.5) for 
incoherent neutron scattering by an ideal gas. This is 
done in the next section. The point to be made here 
is that the decomposition of the memory function 
implied by eq (3.7) is not particularly useful for 
evaluating K{t) , at least for this example. 

4. Ideal Gas Memory Function 

In this section we examine the memory function for 
incoherent neutron scattering by an ideal gas. The 
phase function u is exp(zq-r) where q is the momen- 
tum transfer and r is the position of a gas atom. It is 
a simple matter to show that for the ideal gas 



(cos (q • \t)) = exp [-<7 2 £ 2 /2/3ra] 



(4.1) 



where m is the mass of the gas atom and (3 = l/k B T, 
the inverse of Boltzmann's constant ks times the 
absolute temperature T. 

The memory function is formally 

K(t)=(e H « r (-iq-t;)e-( 1 - p > l ^(iq-v)e l ' qp >. (4.2) 

The evaluation of this expression is now to be under- 
taken. The procedure we follow is to expand the 
exponential factor exp[— i(l — P)Lt], apply the oper- 
ator — i(l — P)L the appropriate number of times and 
then collect the resulting terms in such a way that a 
closed form is obtained for the Laplace transform of 
K(t). It does not seem to be possible to obtain a closed 
form for K(t) itself. 

For the ideal gas, the dynamics are trivial; 



— iLv = v 

- iXv = 0. 

Let us consider first the expansion of K(t) : 



(4.3) 



K(t)=l e-" r (-iq-v) £ 

(-tna-p)iL]" 



n = 

(tq • v)e"«' r 



(4.4) 



We note that the odd terms in t will vanish when the 
expectation value is taken. First consider the n = 2 
coefficient. 



[(l-P)iL] 2 (iq -vW** 
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= (1 -P)iL[(q • v) 2 -<(q • v) 2 )]e i< " r Thus we find 

= [(q • v) 2 - <(q • v) 2 )](-iq • v)e«- K(t) = (x 2 > +C g, (t) 

-M-*i-).'-. <4 ' 5) -((*-<*» J | g Ck <^)^^). 

Repeated application of [(1 — P)iL] 2 yields a sequence j 1 n j+i fc o 

of terms of the type Min 

(4.11) 

( — 1) ten Uq * v)e iq r jyj ow reorc | er t h e y an( j £ sums anc i extract the A; = 

, , „. . . . r term to obtain 

where the coefficients km satisly 

«2n + 2 = A> 2 „(q * v) 2 — (£2/1 (q * v)-) (4.6) 

and J {X 2- (X 2 }) ^ J J C,(^-^)^^^- 2 V 

A: = 1. N fc=lj=fc+l n=j+l \^ n )' I 

A more useful grouping of terms is achieved by in- (4.12) 

troducing another set of coefficients {C,} defined by 

This process can be repeated with the result 

k 2n = 2 Cj(x 2 "-« - (^ 2 "-^>) (4.7) 

K(*) = <* 2 >-Co J ("1)^(0 (4.13) 

iV=l 

with Cq= 1. The notation x= q * v has been introduced w here 
to simplify the writing. The C/s satisfy the recursion 

relation -. **<«)-/<*•-<*■» i i 

C» = -2^'*- 2j >. (4.8) \ *£,* 2 =^ + . 

(— j2)n 

Now let us express K(t) in terms of the C/s. The ... ^ /o ^ t (* 2Al )(x 2k 2~ 2k i) 

factors exp(iq-r) yield 1 when the expectation value n=k N _ t +\ \^ n > - 

is taken so we drop those terms. * 

. . . (x 2k N-l- 2k N-2)(x 2n - 2k N-l\- (4.14) 

\ ^ () (zn) ! / further progress is possible it we make use of the 

moment property of the *'s for the ideal gas, namely 

= <* 2 >+(* 2 i-7^rT£ Cj&—i-(x*—J))) {xV) = (2j-\)\\(x*)K (4.15) 

n=l V ' ' j=0 ' 

,4 g\ The (2n) ! term in the denominator can be removed by 
taking Laplace transforms [15]. For the case N = 1, we 

It will prove convenient to reorder the double sum in 

eq (4.9) so that = ( x 2 - (x 2 ) ) V ( ~ 1)W * 2n \ 

K(t) = (x 2 ) + ((x 2 -(x 2 ))f J Cjt^x 2 "- 2 ) = lx*-x 2 {x 2 ) _\ \ 

/ r 2\ . r A r 2_/ r 2U V ( ~* 2) " ,-A /4,1m Closed forms for the other ftW are also Possible. 

(xj-r-Co^l* U;;2, (2n)! * / (4 * 1U) We evaluate g 2 (z) to illustrate how this is done. 

First we use eqs (4.14,4.15) to obtain 



t 2\n 



+ {<*•-<*-»! I ^l^T- 2 "- 2J ) • &(*)=( (*-<*'» J J (2,-1)!! 



J =1 "=-> +1 j=l n=j+l 

The next step in the analysis involves the use of eq(4.8) (x 2 ) j ( — l) n x 2n ~ 2j \ 

to systematically replace the C/s in eq (4.10) with C () . z 2n + 1 / 
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( * 2 z ( * 2)) t ( 2 >-D !! 



(x 
x 



2 ) "V (-X 2 /2 2 )^ +1 

2 J l + * 2 /z 2 



" 7 *(i /iff) ) | W- D»(-i)^'[WW • 



Next we make use of the result 



= V-J o dq q* exp(-q*l2) 



(4.17) 
(4.18) 



to obtain 



vhere 



:"^tSt» 



iw-^r*?* w 




77 Jo 4 Z 2 +<* 2 >(? 2 



(4.20) 



It is a straightforward task using these techniques to 
show that 

&(*)=(-l)" +1 gi(*)[/(z)]*-». (4.21) 

Next, the Laplace transform of the memory function 
is found to be 



km-— + ww 

This in turn is shown to be 



(4.22) 



£(z)=-z + 



L Vtt<* 2 



Xexp (z 2 /2<% 2 >) erfc (z/V2<* 2 > ) 



(4.23) 



One may also express K(z) in terms of i//(z) as 

K{z)=-z+l$(z) (4.24) 

by taking Laplace transforms of eq (2.3) and solving for 
K(z). In terms of (x 2 ) , it follows from eq (4.1) that 



<M0 = exp {-(x*)t*l2) 



(4.25) 



since 



r.dv{qv) 2 e-* mv *l* 

(x 2 ) =^ = q*lpm (4.26) 



/; 

Jo 



dve-e™ 2 ' 2 



and thus 



*U) 



77 



2<* 2 > 



exp(z 2 /2<* 2 >) erfc (z/V^x 2 )) (4.27) 



in agreement with the result of summing K(z) . 

Although a closed form for K(t) has not been 
found, we have determined the underlying structure 
of the memory function. How it compares with the 
structure implied by eq (3.7) is worth noting. From 
eqs (4.21) and (4.22) it is evident that 

K{z) = ^ + gl (z) f [-ft*)]*. (4.28) 

2 n = 



(4.19) and 



A simple calculation shows that 

gi(t)-c(t)-<**>*(0 

I(t) = d$(t)ldt. 
An iterative solution of eq (3.7) yields 



(4.29) 
(4.30) 



K(z) = C(z) ^ [-/(*)]"■ (4.31) 

Consistency requires that 



fa) £ [-/(*)]» = l/z. 



(4.32) 



Since I (z) — zi//(z) — 1, eq (4.32) is satisfied and the 
two solutions are equivalent. The mathematical 
structure which was found to be useful for summing 
the memory function expansion may appear in retro- 
spect to be inconvenient because of the extensive 
cancellations which necessarily occur. However, 
attempts to impose the structure of eq (3.7) on the 
expansion were unsuccessful in achieving a form which 
could be completely summed. 

In order to see just what sort of structure we have 
developed, let us review the steps followed in summing 
the expansion of K(t). The set of coefficients {kzn} 
resulted from grouping the sum in terms of powers of 
the time. Next, the sum was arranged in terms of the 
deviations of x 21 from the average (x 21 ) and resulted 
in the set of coefficients {Cj}. This was the crucial 
step in the analysis as it exhibited the underlying 
structure of the memory function. Clearly this way of 
arranging the sum is different from the structure 
implied by eq (3.7) since C(t) does not involve the 
projection operator. Next we eliminated all C/s except 
Co from the expansion through repeated use of eq 
(4.8). The closed form obtained in eq (4.23) required 
an explicit relation among the moments of x, namely 
eq (4.15) which follows from the Boltzmann distribution 
for the velocities. It would be interesting to see if this 
type of structure could usefully be imposed on K(t) 
for an interacting system. This topic is not considered 
here. 



5. Discussion 



The time 
figures 1-4 



correlation 
typical 



are 



functions, 
of the on< 



ifj{t), shown in 
;s obtained from 
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liquid state scattering and absorption data. The curves 
initially exhibit a rapid falloff and then go over to 
an approximately exponential function for long times. 
The functions, \p(t), are everywhere positive and, 
except for the short time part, the curves are con- 
cave upward. The memory functions K(t) exhibit a 
rapid initial falloff which carries them to negative 
values. The negative portion of K(t) has the effect 
of decreasing the magnitude of difj/dt and thereby 
reducing the slope of \fj(t) so that ifj(t) remains 
positive for all times. 

It is interesting to compare this situation with that 
found from computer simulations for the velocity 
time correlation functions in Lennard-Jones liquids 
[9, 10]. There the time correlation functions achieve 
negative values. This implies that the initial state is 
"strongly remembered" by the system, a viewpoint 
consistent with the behavior of K(t) since K{t) 
remains positive. At higher temperatures, where the 
kinetic energy is large compared with the strength of 
the pairwise interaction energy, the time correlation 
functions do not go negative while the memory func- 
tions do. 

This suggests that the persistence of the memory 
function is closely related to the time evolution of 
the cage-like many body structure described by 
Rahman [16]. He found from his molecular dynamics 
calculations that a liquid atom undergoes a strong, 
velocity reversing collision and then "rattles" about 
in a loose cage of neighboring atoms until it escapes 
when the cage breaks down. The rattling motion leads 
to negative values of the velocity time correlation 
function. However, if the cage breaks up too rapidly, 
as it would at high temperatures or low densities, no 
rattling motion, coherent with the initial condition, 
occurs andi//(£) remains positive. 

The calculation of the memory function for the ideal 
gas shows that this function has an extensive mathe- 
matical structure. This is expressed in the properties 
of the coefficients {km} and {Cj}. The memory func- 
tions constructed from data show more structure than 
the memory function determined numerically from the 
ideal gas time correlation function. 

This comparison is the basis for the following con- 
jecture. The mathematical structure of the memory 
function for an interacting system should exhibit two 
parts, in contrast with the one part exhibited by the 
ideal gas. There should be a "short-time" part which 
dominates the initial behavior of K(t) and then goes 
rapidly to zero. The second part of the memory func- 
tion should dominate the remainder of the time in- 
terval. This part can be either positive or negative 
depending on the type of correlation involved and is 
dominated by many-body effects which are only poorly 
understood. The testing of this conjecture is a subject 
for further work. 



differential equation by a set of coupled difference 
equations defined on a discrete set of evenly spaced 
time points U with ti +i — ti=kt. For simplicity, a sub- 
script notation is used so iff(ti) becomes i/// and 
i//(0) = i// . The trapezoidal rule and simple difference 
expressions are used to replace integral and differential 
operations. A typical member of this set of difference 
equations is 



1 
2At 



[t//i + i-i/>i-i]=-A* 



1 



(Kofa+Ki^ 






(6.1) 



These equations for K\ are solved recursively in terms 
of the {ijjj} subject to the initial conditions 



Xo = 2(^-i/io)/(A0 2 . 



(6.2) 



The last condition follows from eqs (2.5 and 3.2) which 
indicate that 



(6.3) 



(6.4) 



K(O)=-d*0(*)/d* f |i.o. 

The solution has the form 

Ki = -Koilfi- [i// 2 -i//o]/(A0 2 



i > 1 
-[*i + i-*«-i]/"(Ai)*. 



These expressions were used to obtain the K{t) 
curves discussed in the text. 

It is necessary that \\)(t) be a smooth function of t 
if this technique is to yield useful results for K(t). 
Figures 1-4 indicate the kinds of uncertainties which 
show up in K{t). Some care is needed in the selection 
of the time increment A*. If it is too large the difference 
expression for the first derivative is inaccurate. If Af 
is too small, the noise in \\f(t) governs the outcome 
rather than i//(j) itself. 

Another procedure for solving for K{t) is described 
in appendix B of reference [6]. This procedure is more 
accurate than the one described here in that it uses 
spline fits to ifj(t) to estimate derivatives and it uses 
a more accurate integration scheme. However, when 
the procedure we use is stable, the results of the two 
methods agree closely. 



6. Appendix 



The numerical inversion of eq (2.3) to obtain K(t) 
from \jf(t) is accomplished by replacing the integro- 



I wish to thank Drs. Livingston and Dardy for making 
their data available to me prior to publication. 
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